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Abstract 

We present results for the photon spectrum emitted in non-hnear Compton scattering of pulsed 
ultra-strong laser fields off relativistic electrons for intensities up to ao ^ 100 and pulse lengths 
of a few laser cycles. At ultrahigh laser intensity, it is appropriate to average over the sub- 
structures of the differential photon spectrum. Supplementing this procedure with a stationary 
phase approximation one can evaluate the total emission probability. We find the photon yield in 
pulsed fields to be up to a factor of ten larger than results obtained from a monochromatic wave 
calculation. 
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I. INTRODUCTION 



The ongoing progress in laser technology has opened an avenue towards studying strong- 
field QED effects in ultra-intense laser fields. Presently, the strongest available laser systems 
have a few petawatts, with a focused peak intensity in the order of / = 10^^ W/cm^ pj. In 
the near future, the high-intensity frontier will be pushed forward with ELI [2] exceeding 
/ = 10^^ W/cm^, eventually. A Lorentz- and gauge-invariant dimensionless parameter char- 
acterizing the intensity is given by Og = 7.9 x 10^"'^^/[W/cm^]A^[/im] [3] where /[W/cm^] is 
the laser intensity in W/cm^ and A[/im] denotes the laser wavelength in /xm. Thus, values 
of Oq of several hundreds can be reached in the near future. 



Among other topics, the formation of QED avalanches has been discussed recently |H E] in 
ultra-intense laser fields, where a seed-particle leads to the formation of a cascade by con- 
secutive photon-emission via non-linear Compton scattering and subsequent pair-production 
processes. Here, we focus on the non-linear Compton scattering process in laser fields with 
ao ^ 1 with emphasis on finite-pulse envelope effects. The Compton process has been 
studied for moderately strong laser pulses ag = 0{1) both within classical electrodynamics 
as non-linear Thomson scattering [BHO] as well as in quantum electrodynamics pUl - fTH] . In 
Oq = 0{1) laser pulses, with a duration of several cycles of the carrier wave, the non-linear 
Compton spectrum has interesting structures with many subpeaks per harmonics, which 
may be verified experimentally with present technology [9l [H]. We present here a method 
for calculating energy-averaged photon emission spectra for ao 2> 1 laser pulses. 



Our paper is organized as follows: In section [IT] we evaluate the matrix element for non-linear 



Compton scattering. In section III we use this matrix element to calculate the emission 



probability, presenting our result for the phase-space averaged photon yield. In section 



IV we present our numerical results before concluding in section |V| In appendix |A] we 
comment on the slowly varying envelope approximation. We derive an exact expression 
for the non-linear phase including a carrier envelope phase which has the slowly varying 
envelope approximation an a limit. Finally, in appendix |B] we have collected details on the 
stationary phase approximation, which is at the heart of our approach. 



II. MATRIX ELEMENT FOR NON-LINEAR COMPTON SCATTERING 



A. Basic Relations 



The S matrix element for non-linear Compton scattering — where an electron with mo- 
mentum p emits a single photon with momentum k' and leaves the interaction region with 
momentum p' — in the Furry picture is 



Sfi = —ie j df'x'^pi{x)f'^p{x) 



(1) 



where the Volkov state '^p{x) [15j is a solution of the homogeneous Dirac equation in an 
external classical electromagnetic plane- wave background field with four-potential ^^^(0), 



+ m)*p(x) = 0. 



The Volkov solution reads 

^p{x) = exp{-ip ■ X - ifp{x)} 
where the non-linear phase function is 



1 + 



2k -p 



X) 



Ur. 



(2) 



(3) 



1 



2k -p 



b'[2ep- A{<p')-e^A\(t)')\. 



(4) 



In light-cone coordinates, e.g. x-t = x° ± x^, x_l = and x = (x+,x_l) for the spatio- 

temporal coordinates, the laser four-momentum k has only one component which we choose 
to be k_. We also define a special reference frame (s.r.f.) in which the electron is initially at 
rest, i.e. it has four-momentum p = (m, 0, 0, 0). Assuming a head-on collision of the electrons 
with the laser pulse, the laser frequency in the s.r.f. is u; = (a/7o — 1 + lo)^^o, where uq (70) 
is the laser frequency (Lorentz factor of the electron) in the laboratory system. For the 
vector potential A'^ we employ a transverse plane wave, modified by an envelope function g 
with pulse length parameter r. 



A'^ = AQg{(j)/ t) (e^ cos ^ cos + 63 sin ^ sin > 



(5) 



with linear polarization vectors ej ■ ej = —Sij and ei ■ k = 0, i, j E (1,2), ^ denotes the 
polarization state of the laser [T2j. A complex circular polarization basis e'^ = cos^e^* ± 
z sin 1^62 is suitable for the following considerations. An infinite monochromatic plane wave 



is recovered by the limit (7 — )■ 1. We require that g is a symmetric function of with g{0) = 1 
and g{±oo) = 0. The dimensionless laser amplitude Oq is defined with respect to the peak 
value of the vector potential as = e^o/m. 

The matrix element Q reads in the slowly varying envelope approximation (see [12] for 
details), where terms of 0{t~^) are neglected in the exponent (cf. Appendix [A|) 



S,^(2.)3.'(p'.k'-p)A_l_.^, (6) 
where we have separated the Dirac structures 

^ Ak-pk- p' ' ^">^P'^^P' = . pk . p'^P'^^^^ ^/±'"p' (9) 

from the phase integrals 

00 

^-(s) = [ ci0^'™(0)e'(^-")'^-^^W, (10) 



with /(0) = fp{4>) — fp'{4') and dimensionless momentum transfer 

_ k'_ + p'_ - j9_ s£.f. cy 1 

^~ kl ~ '00 1- ^(i + cos^)^ 



(11) 



where the last equality holds in the special reference frame. The variable s might be in- 
terpreted as a continuous number of absorbed photons, as advocated e.g. in [121 [IS]- The 
reasoning is based on the fact that the energy-momentum conservation can be written com- 
pactly in the suggestive form 

p + sk = p + k'. (12) 

Since the frequency uj' of the perturbatively emitted photon has to be positive, also s has 
to be positive; = would imply s = 0. This holds in any reference frame due to Lorentz 
invariance. 



The non-linear function /p(0) in the exponent in (10) can be split into a rapidly oscillating 



contribution / and a slow ponderomotive part ((/)), i.e. f = f + ((/)), where the latter one is 



averaged over the fast oscillations of the carrier wave. One finds for these two contributions 

/ = ag{<P) sin(0 + 0o) + ^^^g\<P) sm20, (13) 
((/))= /3 [\<P'g'{<P')=:fiGM (14) 

with the coefficients 







a = y al + al = aos sin 6'y cos^ ,^ cos^ + sin ^ sin ip, (15) 

m^an / 1 1 \ srf On/ „N / N 

where the r.h.s. equalities hold in the special reference frame, and (po = arctan2(— 0:2, ai). 
The angles 6 and ip are the polar and azimuthal angle, respectively, of the direction of 
the perturbatively emitted photon, measured w.r.t. the z-axis. Intermediate definitions to 
arrive at the given expressions for a and /3 are ai = ei ■ PcosC,, 02 = €2 ■ Psin,^ and 
= mao {p^,/k ■ p - p'^/k ■ p'). 



B. Evaluation of the oscillating part and harmonics 

In order to calculate the matrix element ^ the integrations in ( 10 ) have to be performed. 



This can be done by a direct numerical integration as demonstrated e.g. in [TT1|T2] for certain 
ranges of parameters. However, if the exponent is very large, i.e. for oq ^ 1, this option is 



inappropriate and one has to employ different methods for evaluating (10). 



The non-periodic oscillating exponential in ( 10 ) can be expanded into a Fourier series over 
the interval [0 — vr, + vr] with the (/)-dependent coefficients [lOj according to 

0+7r 



00 

(p — TT 



J2 B,{^)e-^'t B,{<P) = ^ [ #'e^^^'-/(^'). (17) 



In the general case of arbitrary elliptical polarization, the coefficients are two-variable one- 
parameter Bessel functions [T7] 

00 

5,(0) = Mag,Pg'cos2a2-Ao) = Yl ^^-2s(a^7) J.(/3/ cos 2^ /2)e-^('^-2^)<^o. (18) 

s=—oo 

We focus here on the case of a circularly polarized laser pulse (i.e. cos 2^ = 0), where the 
coefficients Bi simplify to ordinary Bessel functions of the first kind multiplied by a phase 
factor Be{(j)) = e-'^''"'Je{ag{(j))). 



After Fourier expansion, the matrix element ([T]) turns into a sum over partial amplitudes 



E 



(19) 



{s-e) + %^c^,{s + e'-^^^/ctils -£)+ e-'*o^_\Q'+i(^5 - i) (20) 



with purely real coefficients 



d4>J,_^{ag)g{4>rexp{i{s - £)4> - ifSG^m 



(21) 



with G2{4>) = /q d(f)' g'^((f)') from (14). In the expansion (19), the integer i is the net number 



of absorbed laser photons by the electron thus labeling the harmonics. The seeming con- 
tradiction of continuous s and integer i can be resolved easily: In the energy momentum 



conservation (12) the photon four-momentum k is calculated with the central frequency oj. 



However, the pulsed laser field has a finite energy bandwidth Au/u oc 1/r. Therefore, it ap- 
pears as if there would be a continuous photon number. Also the emitted photon spectrum 
is a continuous spectrum. It is not possible to write the energy-momentum conservation as 
p + ik = p' + k' . Consequently, each of the harmonics i has a broad support, not only on a 
delta comb, as for monochromatic bandwidth-free background fields. 

As a consequence of energy conservation, the absorbed photon number has to be positive. 



> 0, since otherwise the energy u' would become negative. Thus, the sum in (19) starts 



at £ = 1. The precise location of the support of the harmonics is determined by the regions 
where stationary phase points (p^ of (21) exist on the real axis, i.e. for real solutions of the 
equation 



9{ 




(22) 



which are always found as pairs ±0^ where we define 0^ > 0. In (22), both the argument of 



the square root has to be positive and the value of the square root has to be larger than zero 
and smaller than unity due to the assumptions made for the pulse envelope g. Independent 
of the explicit shape of the pulse the allowed range of s — £ is bounded by/3<s — £<0 
(note that /3 < 0). Outside this region the stationary phase has an imaginary part leading 
to an exponential suppression of the coefficient functions c^„(s — i). Explicit solutions for 
0^ are listed in Tab. [ijfor the hyperbolic secant pulse g{(p) = 1/ cosh{0/r} and the Gaussian 



hyperbolic secant Gaussian 



9"{0) 

G^U) rtanh^ ^rerf^ 



2 



</>•*• ''"Arcosh-y/ ^ r-y/ln — 



TABLE I. Explicit forms of the pulse shape functions g, their derivatives and the ponderomotive 
integrals G2. Relations for the stationary phase points 0* for these pulse shapes are also provided. 

pulse (7(0) = exp{— 0^/2r^}. In particular, in our numeric calculations we use the hyperbolic 
secant pulse. 

When passing from distant past, = —00, to the distant future, = +00, the coefficients 
pick up a ponderomotive phase shift 

A/ = hm ((/)) - hm ((/)) = hPr (23) 

where the proportionality factor h depends on the explicit shape of the pulse (e.g. /i = 2 for 
the hyperbolic secant pulse and h = y/n for a Gaussian pulse). 



C. Monochromatic limit 



To contrast the realistic case of a finite-duration laser pulse with the idealized case of a 
infinitely long monochromatic wave we consider the reduction of the former to the latter. 



In the limit of a monochromatic plane wave, g ^ 1, the coefficients (21) condense to 



(s-£)n'(27r)5(s-£-/3)J,_„(a) 



(24) 



with support on a delta comb at the locations 



1 + 



^2 



'1 + cos( 



(25) 
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coinciding with the lower boundary of the harmonic support in the case of pulsed laser fields. 
With this, the 5* matrix reads 

% = (27r)^^5^(p + £fc-p'-A:')M^, (26) 

Mt = e-'''^' + %')Ji{a) + ^/e^<^« J,_i(a) + ^_\e-*^" J,+i(a)] , (27) 



2 2 

which coincides with textbook results (e.g. [IS]) and where p = p + -jj^k denotes the 
intensity dependent quasi-momentum of the initial electron with the effective mass ml = 
p ■ p = m^(l + aQ/2). Thus, in a monochromatic wave, the ponderomotive part of the 
exponent provides a constant rate of phase shift ((/)) /3(f) which leads to the build-up of 
quasi-momentum and effective mass. In a pulsed laser field, on the contrary, we find a finite 



total phase shift A/, cf. (23). 



III. PHOTON EMISSION PROBABILITY 

The differential photon emission probability (or photon yield) per laser pulse and per electron 
is defined by (cf. ^) 

dN 1 e'^oj' 



I E ,,jy„,.J -^f- (28) 



doj'dVt 2 QAn^k ■ pk ■ p 

spin,pol 

The cross section for a particular process is obtained by dividing the photon emission prob- 

2 2 

ability by the integrated flux of the laser pulse N^, i.e. da = dN/N^, where = '^'^"^ T^s 
and the effective interaction time T^s is determined uniquely by the pulse envelope 

ns = uj-' / #/(0). (29) 



The interaction time is also related to the total ponderomotive phase shift via A/ = /SuTes. 



In (28), the coherent sum over all harmonics is a double sum 



where we separate the incoherent diagonal contributions from the off-diagonal elements 
which lead to interferences between different harmonics whenever they are overlapping. The 
overlapping of harmonics can happen only for pulsed fields because of the broad support 
of each harmonic; it does not appear for monochromatic laser fields where the support is a 
delta comb. 
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FIG. 1. The differential emission probability for the first three harmonics £ = 1,2 and 3 (from left 
to right panels) in the special reference frame as a function of frequency scaled frequency lxj' /u and 
scattering angle 9. The color code denotes the decadic logarithm of dNi/duj'dQ in units of inverse 
eV. The laser has an intensity parameter of cq = 2 and a pulse length parameter of r = 20. 



One may define the partial differential photon emission probability for a single harmonic by 



du'dVt 



7 



such that 



dN 



2 GAiT^k ■ pk ■ p 

spin,pol 



dNi + interferences. 



(31) 



(32) 



The partial differential emission probability dNi/ doj' dVt is exhibited in Fig. [T] for the first 
three harmonics £ = 1,2 and 3 in the special reference frame as a function of scaled frequency 
w'/w and scattering angle 9. The spectra are shown for = 2 and a pulse length of r = 20 
for a hyperbolic secant pulse shape. In general, the interferences between different harmonics 
are important for differential observables, in particular, for energy resolved spectra, for not 
too high values of ~ 1. There, the substructures of the harmonics, which can be seen 
in Fig. [T| yield interesting spectral information on the non-linear Compton scattering in 
short laser pulses [121 US]- However, for high laser strength, ^ 1, the differential photon 



emission probability (31) is a rapidly oscillating function of the photon energy oj' . The 



number of peaks per harmonic is given by |A/|/7r oc |/3r| and grows, according to (16), 
oc Qq. That means, when measuring the energy spectrum with a spectrometer with finite 
energy resolution, one actually obtains an averaged spectrum. Denoting this average by 



.), we note that 

(dN) = Y^idNe) (33) 



which means that the interference terms in (32) average to zero. Also the total photon 



yield is determined by the diagonal elements alone, e.g. = A^, i.e. the off-diagonal 
elements give no contribution in the case of ao ^ 1. To calculate the photon yield, it 



is necessary to evaluate the coefficients (21). One could solve this problem by a direct 



numerical integration. This works very well for not too large values of and /3. For 
higher laser intensity, e.g. already for ao > 20, /3 can become very large such that the phase 
exponential is in fact a rapidly oscillating function making a direct numerical integration 
very difficult. 

For rapidly oscillating phase integrals, the stationary phase technique can be applied to the 



integrals (21) as done e.g. [TU]. However, even for large there are regions in phase space 



where the stationary phase method is inapplicable: 

1. ) In the vicinity of the non-linear monochromatic resonance at s — i = P, the stationary 

points are located at the center of the pulse and very close to each other. There, the the 
first derivative is almost zero and, therefore, the stationary phase approximation 

tends to diverge. The exponent has to be expanded up to the third order derivative 
of the phase. 

2. ) The stationary phase approximation is appropriate for large ^ 1 only. However, 

in the vicinity of the forward scattering direction 6 = n — i) it behaves as |/3| = 
'^o^T' thus, the stationary phase method can be applied for angles 'd'^ ^ only. 
Therefore, in the phase space regions where is small, a direct numerical evaluation 
may be used. 

Fortunately, these evaluation techniques complement one another, such that they may be 
combined together by suitable matching conditions *Pi^2 (1 and 2 refer to conditions 1.) and 
2.) above) to allow for an accurate calculation of the non-linear Compton scattering spectra. 
The choice of these parameters is motivated in Appendix |b} (We note that this method is not 
restricted to the non-linear Compton scattering process. It can be easily transferred to other 
strong-field processes such as stimulated pair production [19] or one-photon annihilation of 
e'^e'^ pairs [20].) 
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FIG. 2. Differential photon yield dNi/duj'dfl as a function of scaled frequency lo' /u and scattering 
angle for oq = 5, r = 25, and a hyperbolic secant pulse envelope in the special reference frame. 
In the different areas of phase space (I, II and II, which are delineated by dashed curves), different 
methods of evaluation for the coefficients c™ are employed. Parameters are = 2.33 and ^2 = 10. 

We end this section by giving the formula for the phase space averaged photon emission 
probability which reads 

with the averaged coefficients 

2/ > ^1 and \t/3\ > ^2 =1, 

?/ < ^1 and |r/3| > ^2 ^ H, (35) 
\rf3\ <^2 = III 

with y = {s — i — l3)/\j3g"{0)\^. In Appendix |B| we provide further details on the evaluation of 
the stationary phase approximation. The regions I, II and III, where the different evaluation 
schemes for the coefficients apply, are exhibited in Fig. [2j The displayed pattern persists for 
large values of ^ 1. 
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FIG. 3. Partial photon yield in hyperbolic secant pulse with r = 20 (symbols) compared to the 
appropriate yield in a monochromatic wave (thick solid and dashed curves). Matching parameters 
are = 2.33 and ^2 = 20. Left panel: Behavior of the first two harmonics A^i (red diamonds and 
solid curve) and (blue circles and dash-dotted curve) as a function of oq. Right panel: Partial 
emission probabilities Ni for different laser strength ao = 20 (green squares and solid curve) and 
oq = 100 (brown triangles and dashed curve). 



IV. NUMERICAL RESULTS 



For the numerical evaluation, we consider a laser with frequency (in the laboratory system) of 
w = 1.5 eV colliding head-on with an electron beam with energy of 40 MeV. All calculations 
are performed in the special reference frame. The partial photon yield for the first and 
second harmonic per electron in a hyperbolic secant pulse of duration r = 20 are presented in 
Fig. [3} where it is compared to the appropriate probability expected from a monochromatic 
pulse. To get a finite number for the photon yield in a monochromatic field comparable with 
the pulsed laser field we multiply the photon rate Ni in a monochromatic wave with the 



effective interaction time T^q (see eq. (29)), Ni = TegNi, which is a characteristic number 
for each pulse shape. This corresponds to a normalization to the same energy contained in 
the laser pulse. Within this approach, both the pulsed and monochromatic photon yields 
coincide for low oq ^ 1, i.e. in the linear interaction regime, the dependence on the pulse 
shape drops out. For large values of ao, the photon yield in a pulsed laser field is enhanced 
by almost a factor of ten as compared to the photon yield in a monochromatic wave. The 
differences between pulsed and monochromatic yields at large values of ao express non-linear 
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FIG. 4. The total photon yield as a function of oq. The result for pulsed laser fields is shown 
by symbols (red diamonds) being a factor of two above the asymptotic monochromatic result 
(black solid curve, cf. |22j). The black dashed curve shows the sum of the first 140 harmonics in 
a monochromatic wave (compare [22]). The green dotted curve is a rough approximation to the 
asymptotic monochromatic result (see the text). 

finite-size effects. 

In the right panel of Fig. |3| we display the partial photon yields for fixed Oq = 20 and 100 
as a function of the harmonic number ^ including harmonics up to £ = 10^. The cutoff 
harmonic in a monochromatic laser field can be estimated from the behavior of the Bessel 
functions at high index and argument and is given by ^raa.x ~ c^o- While for low harmonics 
the photon yield A^^ is larger in pulsed fields, the ordering of the two curves changes for 
high harmonics, where the photon yield in a pulsed field pulsed fields is smaller than in a 
monochromatic wave. 

Comparing the total photon yields = Yle-^i^ summed over all harmonics up to £ = 10^ 
in Fig. |4| we find that the total photon yield in a pulsed field exceeds the monochromatic 
result by a factor of two for ao = 100. For the monochromatic result we show the sum over 
the first 140 harmonics (dashed curve) and an asymptotic approximation (representing the 
limit flo c>o) being equivalent to a constant crossed field (solid curve, see e.g. [211 122]). 
Additionally, the dotted curve is a rough estimate for the asymptotic probability = 
OiQEDCio^<P [23], where we take A0 = uT^s as the relevant phase interval, and aqED is the 
fine structure constant. For the hyperbolic secant pulse this yields A0 = 2r. 
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V. SUMMARY AND CONCLUSION 



In summary, we calculate the photon yield in non-linear Compton scattering at ultra-high 
intensity ao ^ 1 for pulsed laser fields. We use different methods in different parts of phase 
space, including a stationary phase approximation, where the non-linear phase exponent is 
very large, and a direct numerical integration in phase space regions where this is not the 
case. At ultrahigh intensity, an averaging over the substructures in the energy spectrum 
is included. We find significantly modified partial and total photon yields in pulsed laser 
fields when comparing to the case of monochromatic laser fields. The partial yield for the 
first harmonic can be up to a factor of ten larger than the corresponding monochromatic 
result. For the highest relevant harmonics for a given on the order of £ ~ Oq, the photon 
yield in pulsed fields is typically lower than in monochromatic laser fields. Nonetheless, the 
summed total probability is by a factor of two larger than the monochromatic result which 
is typically approximated by a constant crossed field at large values of ap. This shows that 
the constant crossed field approximation is not as good for pulsed laser fields. The method 
presented here is also be applicable to other processes like stimulated pair production [19] 
and one-photon annihilation fIU\ . 

Finally, we discuss the relevance of radiation reaction. The radiation reaction is relevant if 
the parameter R = oiqedX'^o reaches or exceeds unity, where x = (^ok -plm? is the quantum 
nonlinearity parameter. For the parameters used for our numerical calculations and = 100 
one finds R = 0.04, i.e. radiation reaction effects can be neglected. For even higher values 
of ao, however, radiation reaction effects become important. Furthermore, the fact that 
the quantity related to the emission probability exceeds unity is a hint that multi-photon 
emission, which is related to radiation reaction [23], becomes important at these intensities. 
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Appendix A: Beyond the slowly varying envelope approximation 



In this section, we go beyond the slowly varying envelope approximation. We give a further 
justification of the approximation and comment on carrier envelope phase effects. The 
oscillating part of the non-linear phase exponent has the form 

(Al) 



/((/)) = J dcj) aig cos{(f) + 4>) + Oi2g sin(</) + + /3 cos 2^g^ cos 2{(f) + 
introducing the carrier envelope phase (j). This leads to integrals of the form 



n G (1,2) 



(A2) 



/i = TRe< -(«! -ia2)e*V"2i^i(l,a;a+ l;-x^) 
a 

h = Pt cos 2^ Re j -e^^V"^ 2^^1(2, 2a; 2a + 1; -x^) 



which can be evaluated for a hyperbolic secant pulse with the substitution x = e''^''^ as 

(A3) 
(A4) 

with / = /i + /2, a = (iT + l)/2 and the hypergeometric functions 2-^1- For short pulses 
r^l, a— )'l/2 becomes real. Then, the hypergeometric functions take the limiting values 

' ' ' (A5) 

(A6) 

(AT) 
(A8) 



2^1(1, a; a + 1; -x^) -> 2-^1(1, ^; ^; -x"^) = arctanx. 



2Fi(2, 2a; 2a + 1; -x^) 2^^1(2, 1; 2; -x^ 



l + x2 



such that 



/i = 2r («! cos (j) + a2 sin (j)) arctan exp 
A = t/3 cos 2^ cos 2(f) tanh — . 



r 



Thus, the complete phase reads for r ^ 1 

/((/)) = r(ai cos ^ + a2 sin <^) arctan exp — + r/3(cos^ ^ cos^ (p + sin^ ^ sin^ 0) tanh — , (A9) 

r T 

recovering the result of for a single-cycle laser pulse with linear polarization, ^ = 0, and 

= 0. 

In the opposite limit of long pulses, r 3> 1, a becomes purely imaginary, a — )■ iT/2. The 



slowly varying envelope approximation can be obtained from eqs. (A3) and (A4) by ap- 
proximating 2-^1 (^, na] na + 1; —x^) — )■ 2Fi{n, na; na; —x^) = {x^ + 1)^^ yielding eventually 
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FIG. 5. The phase shift functions Ai and A2 as a function of pulse length r. 



Eq. (13). Thus, our general result eqs. (A3) and (A4) contains both the previous results in 
the slowly varying envelope approximation for r ^ 1 and the single-cycle laser pulses for 
r < 1. 

One of the main differences between the slowly varying envelope approximation and the full 
result is the build-up of an additional phase shift when going from the distant past, = — cxd, 
to distant future, cj) = +00, in addition to the ponderomotive phase shift A/ = 2/3r. In 
the slowly varying envelope approximation this phase shift is equal to zero. In the exact 
expressions for the phase, however, the additional phase shift is non-zero and depends on 
the carrier envelope phase: 

A/i = 2r (^ai cos + ^2 sin 0^ Ai, (AlO) 
A/2 = 4r/3cos2^ COS20A2, (All) 
Ai = jRe(v^(^)-^(^)j, (A12) 



2 1^"' 2 ' ^'2\ 
A2 = ^Re |l + (1 - 2a)m^^) - ^(a)]| , (A13) 

where ip{z) = ^ logr(2;) is the digamma function. Thus, the total phase shift becomes now 
A/ = 2/3r + A/i + A/2, where the first term is the ponderomotive phase shift originating 
from ((/)). 

The functions Ai and A2, depending only on the pulse length r, are depicted in Fig. |5| 
The additional phase is relevant for r < 5 only. For longer pulsed the phase shifts drop to 
zero exponentially fast. This explains why the slowly varying envelope is much better than 
expected; even down to r = 5 although it appears as an expansion in inverse powers of r. It 
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should be clear that this additional phase shift modifies the positions of the stationary phase 
points and, therefore, also the support and number and position of zeros of the coefficient 
functions dr_„. 



Appendix B: Stationary phase, zero convexity approximations and matching condi- 
tions 



1. Stationary phase approximation 



The stationary phase approximation of (21) at a certain value of s is given by the sum of 



the contributions from the two stationary points 



[s - 



exp < i[^s — 



+ j£-nH(-0*))^7™(-0 




exp < —t[s 



2Je_^{ag{<P,))g^{<P, 



27T 



COS 



.s - 



*' 4 



+ z- 



.TT ■ 



(Bl) 



The oscillations of the coefficient functions stem from the cosine term which is a common 
factor in all contributions to the partial matrix element When squaring the matrix 

element to calculate the partial differential photon emission probability dNg^ we obtain a 
cos^ behavior, which averages to 1/2. Because the interferences cancel on average, we may 



replace the cosine in (Bl) by l/\/2 and write for the averaged function 



-t)) = Je^n{(yg{cP,))f 



An 



(B2) 



We now estimate the number of zeros of c™_^ which is also the number of subpeaks in a 
given harmonic. The functions c™_£ are zero when the contributions from the two symmetric 
stationary points interfere destructively; i.e. whenever the phase fulfills the equation 



/3 



4>kg^{4>k 



TC 







(B3) 



defining a series {sk)k=i...K of zeros of c^„(s— £). The number K of zeros is determined by the 
range of values the term in the square brackets can take restricting the largest and smallest 



allowed values of k. Noting that the function in the square brackets in (B3) is a monochro- 



matic dropping function of its argument it suffices to consider its values at ±oo. However, 
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the first term goes to zero when approaching infinity and the second term (including the 
leading (3 term) is ((/)) . Thus, n^k^sLx — ^min) ~ — ((/))(— oo)| and the number of 

zeros K is determined by the total ponderomotive phase shift A/ as K ~ |A/|/7r oc |/3r|/7r. 
This result is a generalization of the findings of [21] within the classical theory of Thomson 
scattering for hyperbolic cosine pulse shapes in the backscattering direction to the quantum 
theory of Compton scattering for arbitrary pulse shapes and arbitrary scattering angles. 



2. Zero convexity approximation 

In the vicinity of = the stationary phase approximation diverges. One must expand 
the exponent up to the third derivative, around the point of zero convexity (i.e. the second 



derivative of the phase in (21) vanishes) yielding the approximation 



ctnis - £) - 27rJ,_„(a) ( -i— ) Ai(-y), (B4) 



'-'-^ (B5) 



|/3^?"(0)H 

where Ai(y) = ^ dtexp{iyt + it^/3} is the Airy function. 

The height and shape of the main peak of the spectrum are determined by the curvature 
of (?(0) at the center of the pulse, where the second derivative vanishes. The main peak 
is the same for hyperbolic secant and Gaussian envelopes, because the second derivative at 
the origin is the same. The only difference comes from the different behavior of the first 
derivative in the stationary phase approximation. This determines the different shape of the 
spectra in the high-energy tail of each harmonic. 



3. Matching Conditions 

Here we motivate our choices for the matching parameters ^1,2- The parameter deter- 
mines the transition between the stationary phase approximation and the zero convexivity 
approximation. We define the matching point hj y = which should be in the range 
1.4 < < 2.3381 because for smaller the two stationary phase points would be too 
close to each other [10] and for larger values the two approximations start to differ too 
much. We use the stationary phase approximation for y > *Pi (region I in Fig. |2]) and zero 
convexity approximation for y < (region II in Fig. [2]). 
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The parameter ^2 determines the matching between the direct numerical evaluation and the 
stationary phase/zero convexity approximation. The approximations for the functions c™„ 
are used whenever |t/3| > and the full numerical results otherwise (region II in Fig. |2]). 
The numerical result was found to be rather insensitive to the explicit value in the range 
10 < *P2 < 50. 
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